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We investigate the scenario where spatially separated parties perform measurements in randomly 
chosen bases on an A''-partite Greenberger-Horne-Zeilinger state. We show that without any align- 
ment of the measurements, the observers will obtain correlations that violate a Bell inequality with 
a probability that rapidly approaches 1 as A'' increases and that this probability is robust against 
noise. We also prove that restricting these randomly chosen measurements to a plane perpendicular 
to a common direction will always generate correlations that violate some Bell inequality. Specifi- 
cally, if each observer chooses their two measurements to be locally orthogonal, then the A'' observers 
will violate one of two Bell inequalities by an amount that increases exponentially with A^. These 
results are also robust against noise and perturbations of each observer's reference direction from 
the common direction. 



I. INTRODUCTION 

Entangled quantum states can yield correlations be- 
tween spatially separated measurements that are incon- 
sistent with any locally causal theory pQ [2] . These non- 
local (or nonclassical) correlations are a resource [3J |3] 
for a range of information processing tasks such as quan- 
tum key distribution HHH], teleportation [51 [TO], certi- 
fication and expansion of randomness [Tl] [T^] and the 
reduction jl3] of communication complexity [14i . 

These nonlocal correlations are correlations of mea- 
surement outcomes. As such, they are not solely a conse- 
quence of entanglement but also depend upon the choice 
of measurements. This point was emphasised by Bell in 
his seminal work [T], where he showed that the perfect 
correlations exhibited by the spin-singlet state do admit 
a simple locally causal model (see also Ref. [Kj and ref- 
erences therein). 

The demonstration of nonlocal correlations typically 
employs carefully chosen measurements whose implemen- 
tation requires the spatially separated observers to share 
a complete reference frame [161 117] . To circumvent this 
requirement, observers who do not initially share a com- 
plete reference frame could share a particular state that is 
invariant under arbitrary rotations of the local reference 
frames [TH]; however, the state preparation involved is 
relatively complicated. Alternatively, they could use cor- 
related quantum systems to establish a shared reference 
frame which can then be used to align measurements [19] ; 
however, this approach is resource-intensive as it requires 
coherently exchanging many entangled quantum systems. 

Recently, it has been shown that such methods are 
not required to demonstrate violations of a Bell inequal- 
ity [50]. In particular, for N spatially separated ob- 
servers that share a Greenberger-Horne-Zeilinger (GHZ) 
state, it was found that most choices of measurements 
lead to nonlocal correlations between measurement out- 
comes [5D]. Therefore distant observers can randomly 
choose measurements that violate some Bell inequality 
with a probability that approaches 1 as iV increases. 
However, the successful detection of nonlocal correlations 



in this scenario requires checking the measurement statis- 
tics against a set of Bell inequalities that grows exponen- 
tially in TV. 

In this paper, we show that for N observers who share 
a GHZ state and are able to perfectly align a single 
measurement direction (but do not share a full reference 
frame), any choice of measurements satisfying a local con- 
straint will generate nonlocal correlations. Furthermore, 
in contrast to the results of Ref. [20] , verifying that these 
correlations are nonlocal only involves testing the mea- 
surement statistics against two Bell inequalities, thereby 
simplifying the verification process. Moreover, we show 
that as N increases, the amount by which the observers 
violate one of the two Bell inequalities increases expo- 
nentially and, in the worst-case scenario, is a constant 
factor below the maximum violation permitted by quan- 
tum mechanics. 

We also investigate the robustness of the above- 
mentioned results and those presented in Ref. [5D] in the 
presence of uncorrelated local noise. We will demonstrate 
that the ability of observers to obtain correlations that 
violate some Bell inequality is robust against depolariz- 
ing or dephasing noise whether or not they share a direc- 
tion. Finally, we show that even if the observers cannot 
perfectly align a measurement direction and so have ran- 
domly perturbed approximations to the common direc- 
tion, they can still always obtain measurement statistics 
that violate one of two Bell inequalities. 

This paper is structured as follows. We begin in SecjlT] 
by illustrating these results in the simplest case, namely, 
when there are two observers who share a Bell state. We 
then outline the generalization to N parties, discuss the 
relevant Bell inequalities and the methods of sampling 
random measurements that are used in this paper. In 



Sec. HI we consider the ideal case in which the N par- 
ties share a state without any noise and can also share 
a reference direction perfectly. In Sec. |IV A[ we relax 
the first assumption and show that the probability of vi- 
olating a Bell inequality is robust against uncorrelated 
depolarizing or dephasing noise. In Sec. |IVB| we also 
show that the probability of violating a Bell inequality 
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is robust against random perturbations in the shared di- 
rection. In Sec. |V] we discuss the imphcations of these 
results and offer some concluding remarks. 



II. PRELIMINARIES 

In this section, we outline the simplest example, 
namely, when two spatially separated parties each per- 
form two binary-outcome measurements on the Bell state 
1$+). For this case, we define the probability of violat- 
ing a Bell inequality and review the results presented in 
Ref. '20] for randomly chosen measurements without any 
shared reference frame. We then show that if both ob- 
servers perform locally orthogonal random measurements 
in the zy-plane of the Bloch sphere, they will always ob- 
tain correlations that violate a Clauser-Horne-Shimony- 
Hoh (CHSH) Bell inequality [1|2T]. We then outline the 
general scenario for N parties and discuss the relevant 
Bell inequalities. We conclude this section by explaining 
the various samplings of measurement bases employed in 
this paper. 



A. A two-party example 

For the two-party case, a verifier prepares many copies 
of the maximally entangled Bell state 
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(|0>|0> + |1)|1)) 



(1) 



where |0), |1) are the computational basis states, and 
distributes one qubit to each of two observers. Both ob- 
servers choose two measurement bases. For each copy of 
the Bell state, the observers each randomly choose and 
perform one of their two measurements on their qubit. 
The observers then send the verifier a list of the mea- 
surement choice (a binary digit Sk) and the correspond- 
ing outcome o^^ = ±1 for each qubit k. The verifier 
uses the lists to determine if the measurement outcomes 
are inconsistent with a locally causal theory. The verifier 
does this by calculating the probabilities p{ol_^ — o^^) (as 
relative frequencies) that the outcomes satisfy oj^ = oj.^ 
given a specific choice of si and S2- They then determine 
the correlation functions 



(2) 



and seek to determine if the correlation functions are 
consistent with a locally causal model. For two parties, 
the correlation functions are inconsistent with a locally 
causal theory if they violate the standard CHSH [21 [21] 
Bell inequality 



However, the choice of the labels for the measurements 
(i.e., which measurement is labelled by Sk = 0) is ar- 
bitrary, as is the labeling of the measurement outcomes 
. Therefore the correlation functions are also incon- 
sistent with any locally causal theory if they violate in- 
equality ([3| after any relabeling of the Sk and/or o'^^ 
and/or the label k. We will follow the terminologies of 
Refs. [22J [23] and refer to two inequalities that can be ob- 
tained from one another through such relabeling as being 
equivalent. 

There are four equivalent inequalities that can be 
obtained from Eq. ^ by mapping (si,S2) to (si,S2), 
(si ® 1, S2), (si, S2 ® 1) or (si ® 1, S2 © 1), namely, 

^ i^^Q^ _^ ^(Q^ ^) _^ ^(^^ _ ^(^^ 1)1 < 2 , 
52^^«"" = \E{0,0) - E{0,1) + E{1,0) + E{1,1)\ < 2, 
Sr'"" = \EiO, 0) + EiO, 1) - Eil, 0) + Eil, 1)1 < 2 , 

^CHSH ^ I _ ^^Q^ 0) ^ ^^Q^ ^) ^ ^^^^ 0) ^ ^^^^ 1)1 < 2 . 

(4) 

All permutations of the o^^ and k map Eq. ([s]) to one 
of the four inequalities in Eq. Q, so the four inequali- 
ties in Eq. Q are the complete set of Bell inequalities 
equivalent to the standard CHSH inequality. This set is 
referred to as the class of CHSH Bell inequalities. For 
two parties, the correlation functions are consistent with 
a locally causal theory if and only if they satisfy all in- 
equalities in the class of CHSH Bell inequalities [21] . 

To see that quantum mechanics predicts violations of 
the CHSH inequalities, first note that for a quantum state 
p and observables O^^ = il^^ ■ a, where a = {(Jx,cry,crz) 
is the vector of Pauli matrices and 



ft': = (sin 6* J cos , sin 6* J sin 0^ , cos 6* J ) 

the correlation functions are 

E{s,,S2) = T^{p{0\®0l^)) . 

For the Bell state p = |$+>($+| defined in Eq. Q, Eq. ([6]) 
becomes 



(5) 



(6) 



E{si,S2) = cos6l,\ cos9% + sin6l,\ sin 6*2^ cos {cj)]^ + <j)%, 



(7) 



If, for example, the measurements correspond to the 
vectors 



r!i = (1,0,0) , 

r!2 = i= (1,1,0) , 



= (0,1,0) , 

n\ = ^{-iXQ), (8) 



then the corresponding correlation functions are 

-71' ^"'"--Tl- 



£^(0,0) 



(9) 



\E{Q, 0) + E{0, 1) + S(l, 0) - E{1, 1)1 < 2 . (3) 



Substituting these into Eq. Q, one finds that 54 
2^/2 > 2, i.e., the CHSH inequality is violated. 
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1. Bell inequality violations with no aligned directions 

From Eqs. Q and ([t]), it is clear that the choices of 
measurements that generate CHSH-inequahty-violating 
correlations must satisfy some constraints, i.e., the di- 
rections that correspond to the measurements must be 
aligned in particular ways. However, if the observers do 
not share a reference frame, directions satisfying such 
constraints can only be chosen probabilistically. If the 
measurement directions $7^^ are sampled according to 
the normalized measure dTl^ , then the observers will 
choose measurements that generate correlations incon- 
sistent with any locally causal theory with probability 

P = y" /cHSH({f^J,f^},f^o,^^?})dr!od^^}dr!2dr!2, (10) 

where /chsh ({^^oi ^^ii ^^o; ^il) ^ function that returns 
1 if the orientations {fij, ilj, J7q, il^} generate correla- 
tions that violate any of the CHSH Bell inequalities and 
otherwise. The probability with which the observers 
generate correlations inconsistent with any locally causal 
theory depends on the way they choose their measure- 
ments, which in turn depends on how much they can 
align their reference frames. For example, if the observers 
can completely align their measurement bases, they can 
simply choose the measurements ([s]), so p = 1, i.e., they 
always generate correlations inconsistent with any locally 
causal theory. 

When the observers cannot align their measurements 
at all and randomly choose both of their measurements 
independently and isotropically, then the probability that 
they will choose measurements that generate correlations 
violating one of the CHSH inequalities is w 28.3% pO] . 
However, if the observers choose their measurements to 
be locally orthogonal, i.e., 

17j-r2j=0, ^1-Vl\ = ^, (11) 

then the probability of generating correlations that vio- 
late a CHSH inequality increases to « 41.3% [20^ 

2. Bell inequality violations with one aligned direction 

Consider another possible scenario, in which the ob- 
servers can align one direction of their measurements, 
e.g., the z direction of the Bloch sphere. Then each ob- 
server can choose two orthogonal measurements in the 
xy-plane, i.e., choose two angles 01 and 02 randomly ac- 
cording to a uniform distribution. The four correspond- 
ing measurements are 

17q — (cos 01, sin 01,0) , 

ill = (—sin 01, cos 01,0) , 

ill = (cos 02, sin 02,0) , 

17} = (-sin 02, cos 02,0) . (12) 



Substituting these into Eq. ([t]) and then Eq. Q gives 

= 2| COS (01 + 02) - sin (01 + 02) I, 
52^«s« = 53°«s« = 0, 

^4™^" = 2| COS (01 + 02) + sin (01 + 02) I . (13) 

Using standard trigonometric identities, we see that the 
CHSH inequalities in Eq. Q are satisfied if and only if 

I cosxl < — ^ , I sina;| < , (14) 
V2 72 

where x — 0i + 02 + f . But one of these inequalities 
is violated unless a; = |^ , so any choice of measurements 
(except for a set of measure zero) will violate one of two 
CHSH inequalities, i.e., p = 1. Therefore in order to 
choose measurements that generate correlations incon- 
sistent with any locally causal theory, it is sufficient to 
perfectly align a single direction, namely the z axis, and 
to check only two Bell inequalities. 

B. The general scenario 

We now generalize the two-party case outlined in the 
previous section to N parties and determine the extent 
to which the probability of generating correlations incon- 
sistent with locally causal theories depends on the align- 
ment of the measurements of the N parties. To this end, 
we consider the scenario (abstracted from the physical 
implementation) wherein a verifier prepares a large num- 
ber of copies of the TV-partite GHZ state (the GHZ state 
is chosen as it is a resource for obtaining maximum vio- 
lations of some commonly used Bell inequalities P51l26| ). 




where |Ojv) and IIat) denote states in which each of the N 
qubits are prepared in the states |0) and |1) respectively. 
The verifier distributes 1 qubit from each copy to iV ob- 
servers. As in the two-party case, each observer chooses 
two measurement bases, which corresponds to the k^^ 
observer choosing two directions il^^ , parametrized as in 
Eq. ([5|, in the Bloch sphere, where G Z2 = {0, 1}. 
Each observer measures their qubits, randomly choosing 
Sk for each qubit. The observers then send a list of the 
measurement labels Sk and outcomes ±1 for each copy 
back to the verifier, who will use the lists to determine if 
the measurement outcomes are inconsistent with a locally 
causal theory. 

In contrast to the typical scenario where {^^s^} are 
fixed a priori to some optimal measurement bases that 
give the maximal violation of a specific labeling of a Bell 
inequality, we now consider a scenario where the mea- 
surement bases/directions {f^s^} are chosen randomly ac- 
cording to some distribution, but are fixed throughout 
the experiment. Formally, if we treat the measurement 
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directions {i^sj.} as random variables, we can define the 
probability, p^, that the verifier identifies that the cor- 
relation functions are incompatible with the class of Bell 
inequalities X as 



where s = X^feLi •^fc and 

P (s, N) = 2— cos {^-{l + N - 2s)j . (20) 

The equivalence class of MABK inequalities is the set 
of inequalities that can be obtained by permutating the 
measurement choices, s^, measurement outcomes, Ok and 
labeling of the observers k in the coefficients /3(s) of in- 
equality (19). However, as we prove in Appendix [B| 



Pi 



N 



k = l SfeSZa 



where d£l ' is the normalized measure associated 
with the sampling of measurement direction 17^'^ , and 
({fig^}) is a function that returns 1 if the measure- 
ments {f^s^} give rise to correlation functions that violate 
any of the Bell inequalities in the class I and otherwise. 

Clearly, depends crucially on both the sampling 
of {f^sfc}j which determines the probability of generating 
nonlocal correlations, and the class of Bell inequalities 
T involved, which determines the probability of the ver- 
ifier identifying nonlocal correlations as Bell-inequality- 
violating. For any given sampling of {fJ^^}, px is thus 
upper bounded by pi^ii=Pnl, where Xaii refers to the 
complete set of Bell inequalities relevant to the scenario 
described above. The feasibility of demonstrating Bell 
inequality violation with randomly chosen measurement 
bases can then be quantified in terms of pnl, which is 
the probability of randomly generating correlations that 
are incompatible with any locally causal theory. We now 
discuss the method of identifying nonlocal correlations 
using an appropriate class of Bell inequalities. 



C. The Bell inequalities 

Bell inequalities are constraints on physically observ- 
able quantities that have to be satisfied by any locally 
causal theory [^. A relevant class of Bell inequali- 
ties for the scenario that we are considering is the set 
of 2^ Mermin-Ardehali-Belinskii-Klyshko (MABK) in- 
equalities [27l [28] . A representative of this class is |29] 



5f =1 E Pi^E{s) 



<2^, 



(17) 



where s = (si, . . . , s^r) is the vector of the N measure- 
ment labels. 



/3(s)= ^ ^/2c 

ae{-l,l}»« 



N 



N 

1=1 



(18) 

a = (ai, . . . , Oat) and the iV-partite correlation functions 
E (s) are the expectation values of the product of the 
measurement outcomes when the fc*'' observer performs 
the Sfc-th measurement. Within quantum theory, the 

N 



maximum possible value of 5^ is 2 5 [25[ [28] . 

As we show in Appendix A, this inequality can be 
rewritten as 



Si 



N 



E 



P{s,N)E{s) 



< 2 



N 



(19) 



all such permuations can be obtained by permuting the 
measurement labels (i.e., Sfc — > 1 — Sk for some set of 
fc e {1, 2, . . . , TV}) and so each of the 2^ MABK inequal- 
ities can be obtained by one of the distinct 2''^ permu- 
tations of measurement settings. In particular, the in- 
equality 



P{N ~ s,N)E{s) 



< 2 



N 



(21) 



which will play an important role in the scenario where 
a single direction is shared, can be obtained from in- 
equality ( 19 ) via the mapping Sk 1 — Sk for all 
fc = 1,2,...,A^. 

When N = 2, the MABK inequalities reduce to the 
Bell-CHSH inequalities [H] and represent the complete 
set of Bell inequalities for this scenario [23]. So for TV = 2, 
the measurement outcomes are inconsistent with any lo- 
cally causal theory if and only if they violate one of the 
MABK inequalities. For N > 2, there are also other 
equivalence classes of Bell inequalities (see, for exam- 
ple, Refs. (25j [SO] [31]). An extensive set of such A^- 
partite Bell inequalities that include the MABK class as a 
subset is the well-known Werner- Wolf-Zukowski-Brukner 
(WWZB) inequalities |25l|29|. These 2^" Bell inequali- 
ties can be put into the form of the following single non- 
linear Bell inequality 



5WWZB 



E 



N 



E 



< 2 



N 



(22) 



Defining Sn = I — N mod 2, the above inequality is 
both necessary and sufficient for the set of 2^ A^-partite 
GHZ correlation functions [with measurements defined 
as m Eq. ^] 



Eis) 



cos 




N 



Sin + Sn II COS e":^, (23) 



k=l 



to be describable within a locally causal theory. However, 
not all measurement statistics are captured by these full 
correlation functions. We can also compute the restricted 
correlation functions of the GHZ state. 



E{{sk}keJc)^Si^l J|cos 0^^, 



(24) 



keK 



/C C {1, 2, . . . , N}, which involve the expectation value of 
the product of the measurement outcomes for a subset of 
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the N parties. As a result, one generally needs to check 
the measurement statistics against the complete set of 
Bell inequalities relevant to the particular experimental 
scenario to determine if these correlations are nonlocal. 
The characterization of the complete set of Bell inequali- 
ties is only known for iV = 2 and 3 (see Refs. [24 l l3U l l3T] 
for details). 

Nevertheless, for small N, the problem of deciding if 
some given measurement statistics are compatible with a 
locally causal description can be solved numerically us- 
ing linear programming For larger values of N, it may 
become infeasible to compute pnl using linear program- 
ming. However, we can make use of the following inclu- 
sion relations: 

{S-f } c {S'f , S^} C MABK C WWZB C Complete Set 

(25) 

to lower bound this probability, i.e., 

P{sf } — -P{5f ,S^} — Pmabk < PWWZB < Pnl, (26) 



orthogonal but otherwise uniform, i.e., RIM with 
the additional constraint: 



where Pmabk etc. are the probabilities defined in Eq. ( 16 ) 
with I being the respective classes of Bell inequalities. 



D. Sampling of measurement directions 

Our sampling of measurement directions depends on 
the extent to which the N observers are able to align their 
measurements within each physical scenario. For exam- 
ple, when all observers share a complete reference frame 
and can completely align their measurements, for any 
class of Bell inequalities I, they can always pick {^^s^} 
in such a way that ({17^^}) = 1, assuming there exist 
such measurements. In this paper, we assume that the 
observers either cannot align their measurements at all 
or can only align them with respect to a single direction 

71. 

The following samplings of measurement directions will 
be applied to cater to the different extents in which the 
observers can align their measurements: 

1. Random isotropic measurements (RIM) - each 
party k chooses both directions fl^^ for = 0, 1 
independently and uniformly from the set of all pos- 
sible directions; 

2. Random orthogonal measurements (ROM) - each 
local pair of measurement directions is chosen to be 



^ The set of locally causal correlations is a convex set with finitely 
many extreme points |32| . To determine if some measurement 
statistics correspond to a member of this set, it suffices to check 
if the given measurement statistics can be written as a convex 
combination of these extreme points. For an alternative, but 
equivalent formulation of this problem as a linear program, see, 
for example, Refs. |33l 134) and the supplementary materials in 
Ref. [20]. 



r^S • f7J = V fc ; 



(27) 



3. Planar random orthogonal measurements (PROM) 
~ in addition to Eq. ( p7| , all measurement direc- 
tions are confined to a plane defined by some nor- 
mal vector n (which corresponds to the common 
direction the observers can align), i.e., 



V fc, Sfe , 
for some n shared by the N parties. 



(28) 



Some of the results presented in Sec. Ill A for RIM and 
ROM have been discussed in Ref. |20j but are included 
here in more detail. 



III. NOISELESS SCENARIOS 

When the N experimenters do not align their measure- 
ments, one may expect that it is unlikely to find Bell- 
inequality-violating correlations by performing measure- 
ments in randomly chosen bases. Nonetheless, for the A^- 
partite GHZ state, the probability of choosing measure- 
ments that violate a Bell inequality rapidly approaches 
1 as A increases. In Sec. |III A| we briefly summarize 
the results for RIM and ROM presented in Ref. |20J and 
analyze the difference between testing the correlations 
against the WWZB inequalities and testing the correla- 
tions against the full set of Bell inequalities (obtained for 
small A^ using linear programming). 

Without any alignment of measurements, if the exper- 
imenters do not test the experimental statistics against 
a class of Bell inequalities that grows exponentially with 
A^, then the probability of identifying that the correla- 
tions generated by the randomly chosen measurements 
are inconsistent with any locally causal theory decreases 
with A. However, if the observers can align the z di- 
rection of their measurements, then we prove that for 
all N they can always choose measurements that vio- 
late one of two Bell inequalities, namely, or S2 from 
Eqns. (19)~(21), by an amount that grows exponentially 



with A^. We also numerically calculate the probability 
of violating four different classes of Bell inequalities and 
show that as the aligned direction is rotated away from 
the z axis, the observers have to test their experimen- 
tal statistics against more Bell inequalities in order to 
certify that the correlations generated by the randomly 
chosen measurements are inconsistent with any locally 
causal theory. 



A. No aligned directions - RIM and ROM 

A natural strategy that the A^ experimenters can adopt 
is to each randomly choose two independent measure- 
ment bases fi^ according to a uniform distribution on 
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0.644% 


15.210% 


25.675% 


99.593% 
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100.0009^ 
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0.428% 


16.879% 


31.235% 


99.965% 




0.734% 


34.442% 


71.190% 


100.000?? 
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0.183% 


19.085% 


37.509% 






0.312% 


41.935% 


80.420% 






10 


0.077% 


20.443% 


42.254% 






0.133% 


47.968% 


86.926% 






15 


0.009% 


22.037% 


50.515% 






0.017% 


59.006% 


95.204% 







There are instances where a randomly generated correlation is local, but our simulation indicates that this happens less than 1 in 
every 10^ times. 

Of the 10^ randomly generated probability distributions, we did not find one that admits a locally causal description. 



TABLE I. Probability of finding a Bell inequality violation from the A'^-partite GHZ state for the scenario where each party 
is allowed to perform binary projective measurements in two randomly chosen measurement bases according to either RIM 
(left) or ROM (right). The number of parties A'' is given in the leftmost column. Then, to the right, we have, respectively, the 

probability of violating Si' , the 2^ MABK inequalities, the 2^ WWZB inequalities, and the complete set of Bell inequalities 
relevant to this scenario. Note that the probability of violation for each class of Bell inequalities is lower bounded by the 
corresponding entry to its left, as expected from Eq. (26 1. 



the surface of a sphere. As can be seen from Eq. ([5]), 
this corresponds to each observer randomly choosing 4 
angles 6*^'^ and (fi'^^ for Sk £ Z2, where <j)^^ are chosen 
from a uniform distribution on the interval [0, 27r] and 
9'^^ from the interval [0, n] with p{9) — ^ sin 6*. When the 
observers restrict their measurements to be orthogonal to 
each other (i.e., when they sample according to ROM), 
then Eq. (27) fixes one of the angles. 



The probabihty of violating 4 classes of BcU inequah- 
ties, namely, {5f }, the 2^ MABK incquahties, the 2^" 
WWZB inequahties and the complete set of Bell inequal- 
ities for two binary-outcome measurements at each site, 
are presented in Tablejlj For N = 2 (and only for N = 2), 
the MABK, WWZB and full set of Bell inequalities are 
all identical to the set of CHSH inequalities, so the prob- 
ability of violating each of these three classes coincide for 
both RIM and ROM. 

For TV > 3, these three classes of inequalities obey 
the strict inclusion relations given in Eq. ( 25 ) , and we 



see that the probabilities of violating these three differ- 
ent classes follow the strict inequalities given in Eq. ( 26 1 



For the MABK and WWZB inequalities, which contain a 
number of inequalities that is exponential in iV, the prob- 
ability of violation px generally increases with N , except 
when N increases from 2 to 3 and a few other cases for 
ROM. This increasing trend is even more pronounced for 
the complete set of Bell inequalities where is found 
to be strictly increasing (up to the limit of our analysis) . 

The WWZB inequalities are necessary and sufficient 
conditions for the full A^-partite correlation functions to 
be consistent with a locally causal theory. Given that the 
restricted correlation functions can be computed from the 
respective reduced density matrices of and are al- 

ways separable, it may seem surprising that the WWZB 



inequalities fail to detect a significant fraction of the 
nonclassical correlations generated from the GHZ states. 
However, while the reduced density matrices of l^'jv) are 
separable and so can be modelled in a locally causal the- 
ory, there is an additional requirement: the locally causal 
models for the different reduced density matrices must be 
consistent, in that they must not contradict one another 
and must also reproduce the full correlation functions 
given in Eq. (23 1. The results given in Table |l] show that 
as N increases, it becomes increasingly difficult to find a 
locally causal model that could simultaneously reproduce 
Eq. (l23l) and Eq. (pll). 



The results from Table |T] also imply that detecting non- 
classical correlations with a probability that increases 
with N requires a set of Bell inequalities containing a 
number of inequalities that is exponential in N . If we 
only use one MABK inequality, e.g., iSf , to detect these 
nonclassical correlations, then the probability of finding 
correlations that violate via ROM decreases exponen- 
tially as N increases. Clearly, because each inequality in 
the same equivalence class can be obtained by adopting 
a different classical labeling, the probability of violat- 
ing any one of the MABK inequalities is equal to 
Therefore the probability of violating one of a set of )vl 
MABK inequalities is upper bounded by MpgN . As pgN 
decreases exponentially with N, M must increase expo- 
nentially with N in order for the probability of violating 
one of a set of M inequalities to cither remain constant 
or increase. 

As we will demonstrate in the next section, this is not 
the case if the N experimenters can align one of their 
measurement directions. In particular, we will show that 
if observers can align a measurement direction, then there 
is a set of two inequalities such that the probability of 
violating either of these two inequalities is one for all N . 
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B. Partially aligned measurements - random 
measurements in the xy plane 

Without any alignment of their measurements, ob- 
servers need to check their experimental statistics against 
an exponentially large class of Bell inequalities to uncover 
nonclassical correlations with a probability that increases 
with the number of observers. However, there are physi- 
cal situations in which it is relatively easy to align a single 
measurement direction, or such a direction is naturally 
defined by the system. 

For example, if qubits are encoded in the polarization 
of single photons and transmitted over optical fibres, then 
the ordinary and extraordinary modes are stable but op- 
tical birefringence causes a phase shift between the two 
modes. If this phase shift is unknown, then the observers 
share a single 'direction' on the Bloch sphere but have an 
essentially random alignment of the other two directions. 
While experimental techniques are available to account 
for this phase shift and may exist for other situations 
in which there is a preferred direction, we show that if 
the observers are trying to violate a Bell inequality, then 
such techniques are unnecessary (the related question for 
quantum key distribution in this situation has also been 
investigated 

Specifically, we show that if the reference direction, 
n, is the z-axis and the observers agree on a labeling 
convention for their measurements, they will always ob- 
tain correlation functions that violate either or S2 if 
the measurements are sampled according to PROM (i.e., 
if the measurements are orthogonal and confined to the 
plane perpendicular to some normal vector n shared by 
the N parties). 

For PROM in the xy plane, the observers share the z 
axis. If the fc*'' observer's two measurements are fig and 
Ui, then, because the labels and 1 are arbitrary, they 
are free to relabel them as necessary so that {n^,n'i,z} 
forms a right-handed coordinate system for all k (a sim- 
ilar result follows for left-handed coordinate systems). 
Under this labeling convention, randomly choosing Hq 
and fli is equivalent to randomly choosing a single ran- 
dom angle Xk, with 6*^^ = | and 



Proof. For the iV-partite GHZ state and the labeling con- 



Xk + Sk - 



(29) 



Theorem III.l. Any choice of orthogonal measure- 
ments in the xy plane on the N-partite GHZ state will 
generate correlation functions that satisfy either 



> 2~ 



„ 3" 1 

> 2^ — ^ 



5f = 2^1 J2 cos((l + 7V-2s)^)£;(s) 
5i^ = 2^| Yl cos[{1-N + 2s)1)e{s) 



(30) 



provided the observers obey the labeling convention de- 
scribed above. 



vention in Eq. ( 29 1 , the full correlation function, Eq. ( 23 1 , 
becomes 

E is) - cos <^ - cos (x + s|) , (31) 

where x = X^fcLi Xk and as before, s — X^^i ^k- Substi- 
tuting this into the left-hand-side of inequality ( 19 ) gives 



5^ = 2"^ I J2 cos((l + 7V-2s)^)cos(x + s|) 



(32) 



There are (^) ways of choosing s such that J2k=i — s, 
so Eq. (l32|) can be rewritten as 



iv+i 1^ /N 



5f = 2^E 



s=0 



(l + A^-2s)^)cos (x + ^l) 



2^^ sin x + (^- 1)- 



(33) 



Similarly, substituting Eq. (31) into the left-hand-side of 
inequality (21) gives 



S. 



N 



Sin (^x + {N + l)j) 



(34) 



Because 



||sin (x-|)|,|sin 



either 5f > 2 



or Si' > 2^-1. 



>_ Vx, (35) 



□ 



From Sec. |II C[ the inequalities in Theorem |III.1| are 
Bell inequalities with an upper bound of 2^ in any locally 
causal theory. Therefore for = 2 and X ¥^ ^ for fc e Z 
the observers will violate or S2 . When N > 2, the 
observers will always violate or 5^ by a factor of at 
least 2^~^. Moreover, the upper bound for both and 
iS^ in quantum mechanics is 2 a ^ so the violation of 
or 5^ is within a constant factor ^ of the maximum 
violation possible in quantum mechanics. 

We can also find the probability p|^max{5;(^, 5^} 

(1 — e)2^T ^ of the observers choosing measurements 

by PROM in the xy plane such that they will obtain 
a violation of a Bell inequality that is within a factor 
(1 — e) of the maximum violation possible in quantum 
mechanics. 

Randomly choosing measurements by PROM in the xy 
plane is equivalent to randomly choosing x in Eqs. (33) 
and (34). The probability of choosing x G [Oj 27r] such 
that 



> 



max{5i^,5f } > (1 - e)2' 



(36) 
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is the same as the probabihty of choosing x € [0, j] such 
that cosx > 1 — e, which is simply - cos~^(l — e). 

Therefore the probability of the observers choosing 
measurements by PROM in the xy plane such that they 
will obtain a violation of either or S2 above some 
threshold (1 — e)2^~^ is 



p(max{5f,5f}>(l-e)2- 



COS 



(37) 



to the set of all normal vectors (i.e., the set of points on 
the surface of a unit sphere with z > 0) by 



Ai = —[\x[ da sin Xg^{a, X) , (39) 
271" Jo Jo 

where 



51 (a, A) 



1 if pi(n) = 1 . 
otherwise . 



(40) 



Partially aligned measurements - random 
measurements in other planes 



Theorem |III.1| applies when the observers measure 
along two orthogonal directions in the xy plane, i.e., when 
the direction that the observers can align is the 2-axis 
(which corresponds to the computational basis used to 
define the GHZ state). When the common direction is 
at some angle to the basis in which the GHZ basis is de- 
fined, the probability of observers obtaining correlation 
functions that violate a Bell inequality can change signif- 
icantly. For iV > 2, we simulate pi{n) as a function of A 
and a where 



(cos a sin A, sin a sin A, cos A) , 



(38) 



and when X is and S2 , the 2 MABK inequalities, 
the 2^ WWZB inequalities or the complete set of Bell 
inequalities for this scenario. 

Given a normal vector, n, shared by N parties, we want 
the probability that n allows the N parties to violate a 
given class of Bell inequalities with probability 1 or with 
some nonzero probability. Consequently, we define the 
ratio of the set of normal vectors that give violations of 
the class X of Bell inequalities with probability px{n) = 1 



180 




270 



270 



FIG. 1. (a) Contour plot of the probability of violating Si" or 
si' for A'^ = 6 when sampling measurements via PROM with 
a reference direction as in Eq. ( 38 1 . (b) Contour plot of the 



probability of violating one of the WWZB inequalities for TV = 
6 when sampling measurements via PROM with the reference 
direction as defined in Eq. ( 38 1. In both contour plots, A = 0° 
in the center and increases radially to a maximum of 90° and 
Pi = 1 for smaller A and generally decreases as A increases 
along a line of fixed a (i.e., along a radial line). 



Similarly, we define the ratio of the set of normal vectors 
that give violations of the class X of Bell inequalities with 
probability px{n) > to the set of all normal vectors by 



Ao = 



1 

2^ 



Jo 



dX / da sin Xgx (a, X) , (41) 



where 



5i(«,A) 



1 ifpi(n)>0, 
otherwise . 



(42) 



Ai gives the fraction of the set of unit normal directions ft 
such that z > and observers who share n can always ob- 
tain correlations that violate a Bell inequality when they 
sample measurements using PROM. Likewise, Ao gives 
the fraction of the set of normal directions n such that 
z > and observers who share n can always obtain corre- 
lations that violate a Bell inequality with nonzero prob- 
ability when they sample measurements using PROM. 
The values of A^ and Ai for N = 2, ... ,6 are given in 
Table [h] and the probability of violating {S^ ,S^} and 
the WWZB inequalities for = 6 is plotted in Fig. [l] 

The case when TV = 2 is exceptional because almost 
any rotation of the reference direction [except when a = 
in Eq. (38l] reduces the probability of violating a Bell in- 
equality to below 1. This occurs because there are dense 
sets of measurements that produce arbitrarily small vio- 
lations of and S2 when ft = (0, 0, 1) and these mea- 
surements do not produce a violation of or S2 (or 
any other Bell inequality) when the reference direction is 
rotated an arbitrarily small amount from the z axis. 

However, for N > 2, there are no such sets of measure- 
ments and so, as our results show, the reference direction 
can be rotated from the z axis by some "small" angle A 
(i.e., A < 35°) in any direction without reducing px. 

Numerically, we found that for any rotated reference 
vector n, considering the full class of MABK inequalities 
provides no advantage over just considering and , 
except when A^ = 2, in which case the MABK inequali- 
ties form the complete set of Bell inequalities. Consider- 
ing the full set of WWZB inequalities does increase the 
value of Ai, but not very substantially. However, testing 
against the full set of WWZB inequalities can substan- 
tially increase the value of Ao, i.e., the area of normal 
vectors for which px{n) > 0, as can be seen in Table |TT] 
Testing against the full set of WWZB inequalities also 
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Ao 




Ai 


N 




MASK 


WWZB 


Complete set 




{5^,52^^} 


MABK 


WWZB 


Complete set 


2 


0.5411 


0.7002 


0.7002 


0.7002 




0.0033 


0.0033 


0.0033 


0.0033 


3 


0.4129 


0.4129 


0.4580 


0.9553 




0.2189 


0.2189 


0.2406 


0.4850 


4 


0.4729 


0.4729 


0.9130 


0.9996 




0.3219 


0.3219 


0.3721 


0.8189 


5 


0.4832 


0.4832 


0.5867 


0.9998 




0.3741 


0.3741 


0.4059 


0.8635 


60 


0.5035 


0.5035 


0.9155 


0.9997 




0.4129 


0.4129 


0.4544 


0.9782 



^ Due to the small sample size for N = 6, we expect the corresponding entries of Aq and Ai to only be, respectively, lower and upper 
bounds. 



TABLE II. Fraction of the surface area of a sphere that corresponds to normal vectors for which the probability of violating a 
class of Bell inequalities is nonzero (Ao) or unity (Ai). The four classes of Bell inequalities considered are {Si' ,82'}, the 2^ 
MABK inequalities, the 2 WWZB inequalities and the complete set of Bell inequalities relevant to this scenario. Note that 
the value of both Ao and Ai exactly coincides for {<Sf ,<S2^} and the MABK inequalities when N ^ 2. 



reveals a strong dependence on the parity of N, which 
occurs due to the 6n term in Eq. (23). 



For N > 2, the dependence on the azimuthal angle, a, 
is small when testing against and 82- In particular, 
P{S" = 1 for all a and A < 35°. That is, for any 

> 3, the reference direction can be rotated from the 
z axis by at least 35° in any direction without affecting 
the probability of generating correlations that violate one 
of two Bell inequalities, namely, and 5^. As N in- 
creases, this threshold value of the polar angle A appears 
to increase slowly. 



IV. NOISY SCENARIOS 

So far, we have made use of various idealisations. We 
now examine what happens when some of these assump- 
tions are relaxed. In Sec. |IV A[ we determine how depo- 
larizing and dephasing noise upon the GHZ state reduce 
the probability of violating a Bell inequality when ob- 
servers do not align any measurement directions or only 
align a single measurement direction. In Sec. |IVB| we 
analyze how the probability of violating a Bell inequality 
is affected by random perturbations in each observer's 
alignment of the common direction, i.e., when the ob- 
servers cannot align their measurements perfectly. 



A. Decoherence 

In order to study the ability of observers to violate a 
Bell inequality in the presence of noise, we consider depo- 
larizing and dephasing noise. For simplicity, we assume 
that each qubit is transmitted over equally noisy, uncor- 
related channels, so the noise for all qubits is described 
by a single parameter v, where v = Q corresponds to 
zero noise and v = 1 corresponds to maximal noise (i.e., 
complete depolarizing or dephasing). We begin by out- 
lining the correlation tensor formalism, which is a con- 
venient method of examining the effect of uncorrelated 
noise. We then give a brief introduction to depolarizing 



and dephasing noise before presenting our results from 
numerical simulations on the probability of violating a 
Bell inequality pi in the presence of noise. 



1. Correlation tensor formalism 

An arbitrary A^-qubit state p can be expanded in any 
basis of Hermitian operators acting on the Hilbert space 

■H2" = (C^)*^^. In particular, the iV-fold tensor prod- 



ucts of local Pauli operators 



(43) 

is one such basis; here a £ Zf'^ is a string of N indices, 
CTaj. G {/2, cr 2^ , CTj^ , CTz } , and I2 is the identity operator act- 
ing on C^. 

Together with the orthogonality relation. 



Tr(EsSj)=2%g, 
we can then represent p by 



(44) 



(45) 



where = Tr [pSj] is the correlation tensor [36 . The 
description in terms of the correlation tensor is thus 
equivalent to the description in terms of the density op- 
erator. In what follows, we will follow Ref. [SS] and de- 
scribe the effect of noise on a quantum state using the 
correlation tensor, which allows us to define the effects 
of uncorrelated noise on each qubit. For the GHZ state, 
we have 



rs-Tr[|«'jv)(*^|Ss] 



= -(OAr|Es|Oiv) 



ljv|S„-|l 



i(6Ar|Ss|ljv) + ^(liv|Sa|0Ar) 



(46) 



All of these terms are unless E3 is a tensor product of 
either (1) 2k Pauli y and {N — 2k) Pauli x matrices or (2) 
2k Pauli z and N — 2k identity matrices for some k G Z. 
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2. Depolarizing noise 
Depolarizing noise maps local Pauli operators as [37] 



(1 - u)ay 



(1 - iy)a^ 



(47) 



Full correlation functions correspond to all observers per- 
forming non-trivial projective measurements, while re- 
stricted correlation functions correspond to some subset 
of observers performing the measurement I2 (i.e., ignor- 
ing the outcomes from some observers). Therefore the 
effects of depolarization on the correlation functions are 



E{s} ^ (1 - iy)^E{s) 



and 



E{{sk}ke.) (1 - J/)l"l£;({sfe}fc6,) , 



(48) 



(49) 



for arbitrary subsets of observers, n C. {1, . . . , N}. 

Note that the effect on the full correlation functions is 
identical to the effect of mixing the GHZ state with the 
maximally mixed state according to 



|*iv) 

when (1 — /i) — 



^ (l-^i)|*^^)(*^^| 

3. Dephasing noise 



T 



(50) 



We also consider dephasing noise, which is appropriate 
when there is some preferred basis in the system which 
is particularly stable. Dephasing noise in the computa- 
tional basis suppresses off-diagonal terms for each qubit, 
i.e., it maps local Pauli operators as [37] 



(1 - v)ay , 



(Jx 



(1 - iy)ax , 

<J, . (51) 



Clearly, from Eqs. (451 and (461 all diagonal terms of 
|^Af)(^Ar| are unaffected and, because off-diagonal terms 
of the correlation tensor are zero unless the term corre- 
sponds to tensor products of only ax and ay matrices, all 
off-diagonal terms are uniformly reduced by a factor of 
(1 — i^)^ . Therefore dephasing takes the GHZ state to 



i^{\0n){0n\ 



iN){lN\)+ ^^ ^' (|OAr)(ljv| + |lAr)(0^|) . 

(52) 

For a dephased GHZ state, the full correlation func- 



tions E{s) oi Eq.(l23l) are 

n 



(l-J.)^COS 




AT 

sin 6'f,+(57v n cos 6*^,, . (53) 



k=\ 



As the restricted correlation functions, Eq. (24 1, depend 



only on the component of the measurements in the z di- 
rection, they are unchanged under dephasing noise (note 
that this is a property specific to the GHZ state). 



Much like its separability property, a decohered or de- 
phased Iv^at) also gradually loses its ability to violate any 
Bell inequality as the level of noise (characterised by v) 
increases. Some bounds on the levels of dephasing and 
depolarising noise at which the I^'at) state no longer vi- 
olates a Bell inequality can be found in Ref. [36] . 



^. Measurements in all directions - RIM and ROM 

For N — 2, . . . , 6 we have numerically calculated the 
probability of violating a Bell inequality, , under the 
influence of depolarizing and dephasing noise for RIM 
(Fig. [2| and ROM (Fig. |3]). For both RIM and ROM, 
the probability of violating a Bell inequality for a given 
level of V is always greater for dephasing noise than it is 
for depolarizing. 

From these plots, we see that for all N , the proba- 
bility of demonstrating Bell-inequality-violating correla- 
tions via either RIM or ROM on the GHZ state is robust 
against depolarizing and dephasing noise. Furthermore, 
this robustness increases with N . For example, with 10% 
dephasing noise (i.e., v = 0.1), the probability of vio- 
lation is reduced to 37% (RIM) or 52% (ROM) of the 
probability in the absence of noise in the case of the bi- 
partite GHZ state whereas for the six-partite GHZ state, 
the probability of violation is reduced to 96% (RIM) of 
the probability in the absence of noise or not affected 
to within the accuracy of the simulations (ROM). Fur- 
thermore, sampling measurements according to ROM not 
only increases the probability of violating a Bell inequal- 
ity in the absence of noise compared to sampling mea- 
surements according to RIM, but also increases the sta- 
bility with respect to both dephasing and depolarizing 
noise. 

It is curious that for some range of noise parameters, 
the bipartite maximally entangled state actually gives 
a higher probability of violation as compared with the 
tripartite GHZ state. 



5. Partially aligned measurements - PROM in the xy plane 

When the observers choose orthogonal measurements 
in the xy plane (i.e., 9 — they will obtain correlation 
functions that violate one of two MABK inequalities by 
an exponential amount. As we now show, this exponen- 
tial violation of a MABK inequality translates directly 
into stability with respect to depolarizng and dephasing 
noise. 

Let us now consider the effect of noise. For mea- 



surements on the xy plane, Eq. (24 1 implies that all re- 
stricted correlation functions vanish. It then follows from 
Eqs. (]48| and (53 1 that the effects of depolarizing and de- 



phasing noise are equivalent and both smoothly reduce 
to (1 — ly)^ . Therefore the observers will always 
violate one of two MABK inequalities with PROM in the 
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FIG. 2. (Color online) Plot of the probability of violation Pnl, sampled using RIM for N = 2, 
of the noise parameter u for (a): depolarizing noise and (b): dephasing noise. 
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FIG. 3. (Color online) Probability of violation Pnl, sampled using ROM for N — 2, 
parameter v for (a): depolarizing noise and (b): dephasing noise. 



, 6 observers, as a function of the noise 



xy plane if 



s AT 3W 

(1 - iy)^2 — - 



> 2 



N 



V2 
V2 



(54) 



Moreover, since the maximum MABK violation attain- 
able by the GHZ state is v2 times the classical upper 
bound correlations generated from the noisy GHZ 
state, Eq. ( [48| , will never violate any MABK inequality 
when 



1/ > 1 



V2 ■ 



(55) 



Therefore, as with the scenario where measurements are 
not restricted to a plane, the ability of observers to always 
violate either or S2 is increasingly robust against 
decoherence as the number of observers increases. For 
iV — >■ 00, these limits are both 1 — -i= and observers will 

(depending on the level of noise ly) either violate or 

'2 



for any choice of measurements or never violate 



or5f. 

The above analysis gives the critical value of v at which 
the probability of violating or S2 reaches the ex- 
tremal values 1 and 0. For intermediate noise levels, 

1 - we can use Eq. ^ 



namely, for G [1 
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FIG. 4. Representation of the normal (n^) and measurement 
directions {^sq and ^sq) for several parties. 



to calculate the probability of violation. From Eqs. (37) 



and ( 48 1 , we have 



p(max{5f (z/)} > (1 - z/)^(l - e)2"^ 

= -cos-i(l-e). 



(56) 



Therefore the observers will obtain statistics that violate 
Si or S2 with probability ^ cos^^(l — e) if 



(1- 1^)^(1-6)2^ =2^ 

^ 1 - e = 



2 — 



N 



(57) 



Substituting Eq. (57) into Eq. (56 1, we find that with 



PROM on the j/-locally-dephased-GHZ state, the ob- 
servers will violate or with probability 



,Sf }(^) 



cos 



2^ 
(1^ 



N 



(58) 



for 1/ e [1 



72 '-^ 



/2 J 



B. Imperfectly aligned measurement directions 

So far, when discussing PROM, we have assumed that 
the observers have a perfectly aligned direction. We now 
consider the case where each observer may have some 
local approximation 

fik ~ (cos afc sin Afc, sinofA; sin Afe, cos Afc) (59) 

to the z axis (i.e., the basis in which the GHZ state is 
defined), see Fig. [4] We assume that the azimuthal angles 



OLk G [0, 27r] are distributed uniformly and the G [0, f ] 
are distributed such that 



P(A') = (1 + T^)[' 



At, 



-^std 



(60) 



Note that with Afc 
convention discussed in Sec 



> ^ are equivalent to with 
< "I, but will change the handedness of the labeling 

for the k^^ 



IIIB 



observer. 



thus changing which Bell inequality will be violated. For 
sufficiently small A, the distribution is analogous to a 
Gaussian distribution with mean and standard devia- 
tion Astd on the surface of a sphere [55] . 

The fc*'* observer then measures in the plane perpen- 
dicular to n/j, i.e., the A;"* observer's two measurements 
are now 



(fi^^ )x — sin cos Afc cos + cos ( 



cos Afc sinoffe — cos 
*: sin Afc . 



smafc 

>t, cos Q!fc 



(61) 



Each observer still chooses random orthogonal measure- 
ments and the same labeling as before, i.e.. 



4^ 



Xk+ Sk- 



(62) 



In Fig. [5] we present numerical results for the proba- 
bility of violating one of the 2^ MABK inequalities as a 
function of the standard deviation of the distribution of 
polar angles Astd for = 2, . . . , 6 observers and the dif- 
ference that checking the full set of 2^ MABK inequal- 
ities rather than just and . Even for standard 
deviations Astd ~ 90°, the probability of violating an 
MABK inequality is higher than when the observers do 
not share a direction, cf. Table |l| The effect of increas- 
ing Astd: that is, of decreasing the average accuracy of 
the local approximation to the common direction, is to 
smoothly decrease the probability of violating a MABK 
inequality. We find that checking against more inequali- 
ties makes little difference for Astd ^ 30° (i.e., the prob- 
ability increases by less than 0.1%). As the individual 
approximations become less accurate, the difference in- 
creases approximately linearly and also increases with N. 
At first glance, this seems rather counter-intuitive, as for 
Astd ~ 90°, the observers have essentially no idea what 
the reference direction is. The difference arises because 
the pseudo-Gaussian distribution is biased to smaller val- 
ues of A, which is only defined modulo 90°. 



V. DISCUSSION AND CONCLUSION 

In this paper, we have shown that the degree to which 
N observers can align their measurements substantially 
affects the probability of them generating Bell-inequality- 
violating (BIV) correlations by performing randomly cho- 
sen measurements on the A^-partite GHZ state. Further- 
more, the difficulty involved in verifying that some cor- 
relations are BIV is also closely related to the extent to 
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FIG. 5. (Color online) (a) Probability of violation pmabk when measurements are sampled using orthogonal measurements 
in the plane perpendicular to the normals Uk that are distributed according to a pseudo-Gaussian distribution as a function 
of the standard deviation of the distribution of polar angles A^td for A'^ = 2, . . . , 6 observers, (b) Increase in the probability of 
observers finding that their correlation functions do not correspond to a locally causal model when they check against the the 
full set of 2^ MABK inequalities, rather than just Si' and 51^ 



which the observers can align their measurements: the 
better alignment they have, the fewer inequalities are 
needed to verify the nonlocal nature of these correla- 
tions. However, this reduction in difficulty does depend 
on how the observers can align their measurements. If, 
for example, they align the measurement directions cor- 
responding to the basis in which the shared GHZ state 
is defined, then the observers will always generate BIV- 
correlations and this can be verified by testing the corre- 
lation functions against just two fixed Bell inequalities. 
As the aligned direction is rotated away from the z axis 
(which corresponds to the basis in which the GHZ state 
is defined), the probability of violating a Bell inequality 
smoothly decreases. These results may make it easier 
to test Bell inequality violations over large distances, as 
they reduce the need to align distant measurements. 

We have also shown that these results and the results 
presented in Ref. [5D] are strongly robust against uncor- 
related noise. Moreover, we have shown that even if the 
observers can only partially align a measurement direc- 
tion, i.e., if each observer has an approximation to the z 
axis that is distributed with a standard deviation of up to 
35° , they can still obtain BlV-correlations with probabil- 
ity 1. This suggests that the idea of demonstrating BIV 
correlations using randomly chosen measurement bases 
is not only an idealization, but is also applicable to real- 
world scenarios. Our results may also have implications 
for reference-frame-independent quantum key distribu- 
tion [35j in the presence of noise. Furthermore, given the 
close connection between Bell inequality violation and 
the security of quantum key distribution protocols, it 



will be interesting to investigate if the results presented 
here have any implications on real-life implementation of 
device- independent quantum key distributions [32! ■ 

Our results also provide insight into the structure of 
the set of locally causal correlations and its relations with 
respect to the set of quantum correlations. To this end, 
consider the results presented in Table |T] As discussed 
in Sec. HI A the probability of violating any one of the 
MABK inequalities is equal. If a given set of correlation 
functions could violate at most one MABK inequality, 
then we would have 



D , Pmabk ,r 
R := log2 = N , 



(63) 



as 5f is one of the 2^ MABK inequalities. For ^ = 2, 
this holds, and one can indeed show that at most one 
CHSH inequality can be violated (see the supplementary 
material of Ref. ^3). However, for larger values of N, 
we see that Eq. (63 1 no longer holds: for example, with 
RIM, we have R = 2.913 and R = 11.258, respectively, 
for = 3 and A'' = 15; similar scaling is also observed 
for ROM. In relation to this, it will also be interesting 
to determine if there arc other aspects about the sets 
of correlations that we can learn by surveying randomly 
generated correlations, a problem that we shall leave for 
future research. 
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Appendix A: Reformulation of the MABK 
inequalities 

In this appendix, we show that the coefficients 



ae{-l,l}®" 



N 



N 

1=1 



'I ' 



(Al) 



in the MABK inequalities depend only on s = X^jLi ■^fe 
and N. In particular, we prove that 

/3{s)=:l3 (s, N) = 2^ cos (^(1 + - 2s)) , (A2) 

for all s e Zf"^ and N > 2. The proof is by induction in 
N. 

Proof. For the purpose of the proof, we will use s^~^ and 
to denote the {N — l)-bit string (si, . . . , sjy-i) and 
the A^-bit string (si, . . . , sn)', likewise for d^~^ and a^. 
Moreover, let us define s' := X^a^/ -^fei 'I2k=i^k, 
Ylk=i o-k and 



7(s 



N 

^ y2sin(^(A + l-a))n«: 



ae{-l,l}®" 



4' 



(A3) 

For N = 2,3, it can be easily checked by explicit calcu- 
lation that Eq. ( A2 1 holds. To prove Eq. ( A2 1 for general 



A' > 2, let us first establish some recursion relations. To 
this end, let us expand Eq. ( A2 1 in terms of ojv = ±1 to 
get 



N-l 



p{s^)= y: ^cos(^(^-a')) n«r 



ie{-i,i}« 



-i-iy- Y ^cos(^(7V-l-2-a')) 

a"-ie{-i,i}®"-i 1=1 

(A4) 



(A5) 



As cos(a: + f ) = ^ sin(a;), we have 

/3(s^)=/3(s^-i) -(-1)^7(5-'^-') • 
Using a similar argument, we have 

7 (s^-\ 1 - ^ 7 G?^-^) + (-1)^-^/3 is''-') 

= -{-ly-p (s^) , (A6) 

or equivalently, 



7(s^)=-(-l)i--/?( 



1 - s 



N 



(A7) 



By the induction hypothesis, let us suppose that Eq. ( Al ) 
is true for some N = no — 1^ where ^ 3, we will now 
prove that it also holds for N — hq. Explicitly, note from 
Eq. (fASl) that 



= /3 (s""-^) - (-1)^-0 +-"0-1/3 (s"o-2^ 1 _ 

= 2^ COS ^(no-2s') - (-1)^-0 +-"o-i2^ x 
(no - 2(s' + 1 - 2s„„_i)) 



COS 



2^ 



= 2- 



4 
cos 

COS 
COS 
COS 



{no -2s') - (-1)^-0 X 



-(no -2s' -2) 



-(no - 2s') 



= 2 



^+1 



cos 



„ "0 + l 

— 2 2 cos 



(no - 2s' + 2 - 4s„J 



(no + 1 - 2s) 
(no + 1 - 2s) 



cos 



"(2Sn|-, If 



where the second equality follows from Eq. (A7|, the 



third equality follows from the induction hypothesis, and 
the other equalities follow from simple algebraic calcula- 
tion using trigonometric identities. □ 



Appendix B: Equivalent MABK Inequalities 

Bell inequalities in the same equivalence class are those 
that can be obtained from one another by some per- 
mutation of the labels on the parties fc, and/or set- 
tings Sfe = o Sfe = 1 and/or outcomes "-|-1" <->■ "- 
1" in the coefficients defining the inequality [52 (see 
also Ref. [25]). Testing a given set of correlation func- 
tions against a different, but equivalent Bell inequality 
amounts to testing the same Bell inequality, but with 
a different labeling and/or sign on the correlation func- 
tions. 

Given that there are A^! permutations on the label fc, 
2^ distinct permutations on the labels Sfc, and two differ- 
ent labeling of outcomes for each of the 2A'^ measurement 
directions, the number of inequalities that are equivalent 
to Eq. (20) is upper bounded by A!2^^. However, as we 



will show, most of these relabelings give identical inequal- 
ities. To this end, let us start by proving the following 
lemma. 



Lemma B.l. For the MABK inequality in Eq. ( 19 1, any 



relabeling of the measurement outcomes on a subset of 
parties k d Af can be .simulated by a permutation of the 
label Sfe for the remaining parties. 

Proof. Firstly, note that the effect of relabeling the mea- 
surement outcomes for the j-th measurement setting of 



15 



the fc*^ party modifies the correlation functions E{s) by 
a phase factor, i.e., 



(Bl) 



From this, we can see that relabehng the outcomes of 
both Sk — and 1 will only introduce a global sign 
change, which has no effect because of the absolute value 
function. Henceforth, we therefore only consider the case 
where the outcomes of the measurement corresponding to 
Sfc = 1 are relabelled. 

For simplicity, let us also consider the scenario where 
only the label for the measurement outcome of the A:*'* 
party is changed (the proof for the more general scenario 
proceeds analogously). We can absorb the effect of this 
relabehng into the /3's by setting 



l3\s,N) = {-iy''(3{s,N). 



(B2) 



We now show that this sign change can be simulated by a 
change of the label Sk for the remaining — 1 observers, 
i.e., by setting 



s[ = Ski s; + (1 - - si) 



(B3) 



With this new labeling, we have 



N 



N 



a-si)^N-l-s + 2sk. (B4) 



Substituting this into Eq. ( 20 ) gives 



Pis\N) = 2— cos i^-{-N + 3 + 2s~4sk))) 
= ^{-IY''2^ COS (^(-^ - 1 + 2s)) 

jV + 1 / TT 

= -(-1)"'=2— cos i^-{N +1- 2s) 



-(-l)^^/?(s,iV), 



(B5) 



which is the same as Eq. ( |B2[ ), except for an overall sign 
which has no effect because of the absolute value func- 
tion. □ 

It is then a small step away to prove the following 
Theorem. 

Theorem B.2. There are 2^ Bell inequalities that are 



equivalent to Eq. ( 19 ) under relahelings of the measure- 



ment outcomes, measurement labels or permutations of 
observers. 

Proof. There are 2^ independent permutations on the 
measurement settings, mapping Sfc = O = 1, which 
give different inequalities since all the s^s are indepen- 
dent. By Lemma |B.1[ we know that the permutation 
of measurement outcomes do not introduce any new in- 
equality beyond this set of 2^ inequalities. 

Next, note that a permutation of observers rearranges 
the label A;, which corresponds to a permutation of the 
entries in the vector s — (si, . . . , sn)- Obviously, this 
leaves the sum of the entries in s unchanged. Because 
j3{s,N) only depends on s and iV, /3 is invariant under 
permutations of the entries in s, so permutations of the 
observers simply rearrange the terms in Eq. (19). □ 
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